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Abstract. We investigate a chain of superconducting stripline resonators, each 
interacting with a transmon qubit, that are capacitively coupled in a row. We show that 
the dynamics of this system can be described by a Bose-Hubbard Hamiltonian with 
attractive interactions for polaritons, superpositions of photons and qubit excitations. 
This setup we envisage constitutes one of the first platforms where all technological 
components that are needed to experimentally study chains of strongly interacting 
polaritons have already been realized. By driving the first stripline resonator with 
a microwave source and detecting the output field of the last stripline resonator one 
can spectroscopically probe properties of the system in the driven dissipative regime. 
We calculate the stationary polariton density and density-density correlations ^^^^ for 
the last cavity which can be measured via the output field. Our results display a 
transition from a coherent to a quantum field as the ratio of on site interactions to 
driving strength is increased. 
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1. Introduction 

In recent years, the investigation of condensed matter and quantum many-body 
systems with quantum simulators, artificial quantum many-body systems that offer 
unprecedented controllability and measurement access in the laboratory, has become 
an active research area and is currently receiving increasing attention. The technology 
employed for quantum simulators ranges from ultra-cold atoms [Ij to ion traps [2j and 
systems of coupled cavities [31 HJ [5l f6j among others. 

Due to recent technological progress, arrays of coupled cavities and optical 
nano-fibers in which the trapped light modes couple to atoms are now becoming 
suitable devices for the generation of quantum many-body systems of polaritons 
[SI Hi |5l El |7[ |8[ [9l [ID], i.e. quantum-mechanical superpositions of atomic and photonic 
excitations. In these systems, it is of particular interest but also most challenging 
to reach a strongly correlated regime, where their dynamics diflFers most significantly 
from that of classical light fields. The key experimental requirement for reaching these 
conditions is a so called strong coupling regime for the cavities respectively the fiber. 
This means that the coherent coupling between the light modes and the atoms or other 
optical emitters must be strong compared to the loss processes which are inevitably 
present in every device. A very impressive strong coupling regime has recently been 
realized in circuit cavities, making these devices an ideal platform for studying strongly 
correlated polaritons. 

Circuit QED [131 El US] was developed as a solid state equivalent to optical 
cavity QED, coupling Josephson qubits that are acting as artificial atoms with stripline 
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resonators acting as cavities for microwave photons. Here, the reduced quasi one- 
dimensional mode volume of the stripline resonator and the enhanced dipole moment 
of the Josephson qubit with respect to atoms give rise to a pronounced strong coupling 
regime, where the coupling between resonator and qubit, significantly exceeds both 
the decay rate of the resonator, and the qubit, g / n ^ 1 and 5^/7 ^ 1. Moreover, 
since stripline resonators trap microwave photons, they are more than 1cm long. The 
precision of current fabrication techniques thus allows to build several resonators that 
can resonantly couple to each other via mutual photon tunneling on the same chip 
[m [m [121 . The currently employed circuit QED technology thus permits to build 
arrays of resonantly coupled cavities that each interact in a strong coupling regime with 
qubit s. In this way it is one of the first setups to feature all properties which are needed 
to generate strongly correlated many-body systems of polaritons. 

In this work we show that an effective Bose-Hubbard Hamiltonian for polaritons 
can be engineered in an array of stripline resonators that each couple to a transmon 
qubit [iTj . Josephson qubits fl^ |T9] come in basically three different fiavours depending 
on the property that is controlled from the outside or rather the channel of the qubit 
environment coupling: fiux-f20]. phase- [21] and charge- qubits [22l[23j. Here we consider 
a setup with transmon qubits [17] which are charge qubits (Cooper pair boxes) operated 
at sufficiently enhanced values for the ratio of Josephson energy, over charging 
energy, Ec^, Ej/ Ec > 50 and are robust against decoherence caused by fiuctuations of 
background charges. We emphasize that the effective Bose-Hubbard Hamiltonian we 
derive can be realized with resonators and qubits of readily existing technology. 

Quantum phases for the ground state and low temperature thermal states of the 
Bose-Hubbard Hamiltonian have been studied with ultra-cold atoms trapped in optical 
lattices [Ij. This system has also been employed to study the dynamics of none- 
equilibrium states that were prepared by sudden quenches of some lattice parameters 
[24] . In contrast, a realisation in an array of stripline resonators allows to investigate 
the Bose-Hubbard Hamiltonian in a fundamentally different regime, where the resonator 
array is permanently driven by lasers to load it with photons and thus compensate for the 
excitations that are lost due to qubit relaxation and cavity decay. Whereas substantial 
understanding of equilibrium quantum phase transitions has been achieved, a lot less 
is known about these non-equilibrium scenarios where the dynamical balance between 
loading and loss mechanisms leads to stationary states. It is the investigation of these 
stationary states, that our approach to the Bose-Hubbard Hamiltonian is ideally suited 
for. 

Experiments with transmon qubits [TTl [25] coupled to a stripline resonator are often 
conducted without directly measuring the state of the qubit but by spectroscopically 
probing the transmission properties of the resonator. In an experiment the effective 
Bose-Hubbard Hamiltonian will thus be operated out of thermal equilibrium in a driven 
dissipative regime [26l [271 [2HI [29]. In a suitable setup with a linear chain of resonators 
one would thus drive the first resonator with a coherent microwave input and measure 
the properties of the output signal at the opposite end of the chain. 
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Figure 1. Sketch of the proposed system to simulate Bose-Hubbard physics. Striphne 
resonators are coupled capacitively in a chain and each stripline resonator is coupled 
to a transmon qubit that gives rise to an on-site interaction for the polaritons. For 
definitions of Jo and see ^ and Q. 



In the regime we consider, this situation can be accurately described by a Bose- 
Hubbard Hamiltonian for polaritons with a coherent driving term at the first site and 
Markovian losses of polaritons due to cavity decay and qubit relaxation at all sites of 
the chain. In this scenario, the interplay of coherent drive and polariton loss leads to the 
emergence of steady states, for which we derive the particle statistics and characteristic 
correlations. In doing so we focus on the polariton statistics, in particular the density 
and density-density correlations, in the last resonator as these can be measured via the 
output signal. Our results show a transition from a coherent field to a field with strongly 
non-classical particle statistics as the ratio of on site interactions to driving strength is 
increased. 

This work is devided into two main parts. In section [2] we show that the dynamics 
of our system can be described by a Bose-Hubbard model and in section [3] we present 
the results of our calculations for the polariton density and density-density correlations. 

2. Transmon-QED and the Bose-Hubbard Model 

To generate a Bose-Hubbard model with polaritons we consider an array of capacitively 
coupled stripline resonators with each resonator coupled to a transmon qubit, see 
figure [ij In this section, we first introduce the Hamiltonian that describes this setup and 
then show how it can be considerably simplified and transformed into a Bose-Hubbard 
Hamiltonian for two polariton species. 

2.1. The full Hamiltonian 

The full Hamiltonian of our setup is a sum of single-site Hamiltonians, i7i_site,i, that 
each describe a transmon qubit coupled to a stripline resonator and terms that describe 
the capacitive coupling between neighbouring stripline resonators, Hcj^i^i+i-> 



A transmon qubit can be regarded as a Cooper pair box that is operated at a large ratio 
of Ej/ Ec ^ 1. This regime can be accessed by shunting the Josephson junction with an 




(1) 
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additional large capacitance and thereby lowering the charging energy Ec = /{2C^). 
Here, Cy, = Cj -\- Cg -\- Cb is the sum of the junction's capacitance, Cj, the mutual 
capacitance with stripline resonator, (7^, and the shunting capacitance Cb- The 
Hamiltonian for one stripline resonator coupled to a transmon qubit reads, 

^i-site = ^Ec (n - nf - n^;Y " ^os ((^) + cj^a+a . (2) 

Here, is the resonance frequency of the isolated resonator and we have omitted the 
site-index i for readability. Transmon qubits consist of two superconducting islands 
connected by Josephson junctions, h is the operator for the difference in the number of 
Cooper pairs on the two superconducting islands, rig^ the offset charge induced by an 
applied dc voltage and intrinsical defects and is the operator for the superconducting- 
phase difference between the two islands. We assume the transmon qubit to be placed 
in the antinode of the stripline resonator's field mode. This gives rise to an additional 
ac component in the offset charge, rig^ = ^9 /(2e)Vjnl{ci + cl^) ^ where K-ms = \fjjffC'r is the 
root mean square voltage of the vacuum field mode and a the annihilation operator of 
photons in the stripline resonator. The offset charge rf^ thus induces a coupling between 
the transmon and photons in the resonator. For circuit QED setups one normally uses 
A-resonators with the antinode located at the middle of the resonator. 

The energy of the coupling capacitor between neighbouring transmission line 
resonators, e.g. sites i and i + 1, can be expressed in terms of the difference in the 
electrostatic potentials across the capacitor, 

_ Cjiyi — ViJ^x)^ _^J^'r ( \ y t \2 

Here, Cr is the capacitance of the whole stripline resonator with respect to the ground 
plane and Cj the capacitance of the capacitor that connects the two resonators. We 
assume the electrostatic potential in resonator i to have antinodes at the ends of the 
resonator and write it in terms of the creation and annihilation operators, aj and a^. 
We now turn to simplify the Hamiltonian ([T]) by a sequence of approximations. 

2.2. Approximations to single-site terms 

We first simplify the single-site terms, i7i_site, as in ([2]). For large Ej/Ec and low 
energies, the phase difference between the two islands remains small and we can expand 
the cosine in ^ around 99 = up to quart ic order. 



^i-lite = 4^c {n - nf - n^;f - Ej + - + uj'^a^a . (3) 



Higher order terms can be neglected, c.f. fTT]. In terms of bosonic creation and 
annihilation operators for the transmon qubit excitations. 



n = 

c 



2 \2E, 
f2Ec 
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the Hamiltonian ^ reads, 

"^1— site -f^transmon H" -f^lin H" -f^coupling H" Hj-^g , (4) 



where 



E 



2^, 



with, 



+ 2^eK(°) (a + at) 
-f^coupiing = ig{h - })^){a + a)) 



= V 8EcEj 



Cj^Cr 

g = %V^'^ ( ^ 



The terms hnear in the creation and annihilation operators can be ehminated by 
performing the unitary transformation 

U =Ui®U2 

Ui exp (ra — ra^) 

U2 = exp {i{sb^ + sb)) . 

that displaces the creation and annihilation operators by the constants r and s 
respectively 

a Uialll = a + r 
b U2bUl = b-is. 

r and s can now be chosen such that all terms linear in a and b cancel in the transformed 
Hamiltonian. Finally the interaction between the transmon qubit and the field mode 
of the stripline resonator gets reduced to an exchange interaction in a rotating wave 
approximation. To justify this rotating wave approximation we have to ensure that the 
interaction strength between the transmon qubit and the stripline resonator is smaller 
than the sum of the frequencies of the two, 

— ^ « 1 . (5) 

Parameters extracted from [30] are cj^ = 43.6Ghz, Ec = 0.4Ghz and a maximal value 
for Ej/Ec of 150. We choose CgcVjnl/iCj^uJr) = 0.1 which is in agreement with the 
theoretical upper bound in [17j and find ^, 0.1. The single-site Hamiltonian can 

thus be approximated by, 

^f-site = V^EcEjb^b -^[b + b^Y + g (at 6 + ab^) + u^a^a . (6) 
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2.3. Approximations to couplings between resonators 

We now turn to simplify the couphngs between neighbouring resonators, Hcj^i^i+i- We 
assume that Cj <C Cr which imphes that CjUj'^/{2Cr) is smaU compared to the isolated 
cavity frequency uj'^^ Cj/{2Cr) <C 1 and apply a rotating wave approximation to neglect 
those terms in the intercavity interaction that don't conserve the total photon number, 

where Jq = §^^r- The first term on the rhs of ^ can be absorbed into the single-site 
Hamiltonians by introducing a shifted resonator frequency 

^r=u';(l + ^+2^(l + p-\] , (8) 

and the remaining term in ([T]) describes tunneling of photons between neigh- 
bouring resonators. Next, we explain how the simplified Hamiltonian H^'^^ = 
J2i (^^i-site i + i i+i) ^c transformed to a two component Bose-Hubbard Hamil- 
tonian. 



2.4' The polariton modes 

In the case of circuit QED with transmon qubits the coupling constant between photonic 
and qubit excitations is the dominating interaction energy of the system. Excitations of 
the whole system therefore can't be characterized as purely photonic or qubit excitations 
in general. To obtain a more suitable description we introduce new creation and 
annihilation operators, 

c+ = cos(0)a + 8m{9)b (9a) 
c_ = sin(e)a-cos(e)6, (96) 

describing excitations commonly termed polaritons where, 

9 



Au + ^Au^ + g^ 
g^+(^Auj + ^Auj^ + g^y 



cos(e) 



with Auj = ujj. — ujq. The sine and cosine terms account for the transition of the 
character of the exitations from photonic to qubit excitations for the c+-mode as the 
ratio Ej/Ec increases and vice versa for the c_-mode. Expressing the Hamiltonian 
^(2) ^ + polariton modes (9a-6) we get, 

7/(2) = i7c+,lin + ^c_,lin + ^cc + ^nlin • 

This Hamiltonian consists of two harmonic chains for the c+ and c_ polariton modes. 
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with (jo'^ = {{ujr +(^q) ± -\/A^~+~^)/2, a term describing hopping from a c_-mode at site 
i to a c+-mode at site i + 1 and all other possible combinations, 

Hcc = - Jo sm{e) cos{0) (^c| +Ci+i _ + c| _Ci+i,+ + h.c.^ , 



and a term describing the nonlinearity, 

-E, 
12 



-f^nlin — 



J2 {MO) {c^,+ + 4+) - COS(^) + . 



We assume the frequencies of the two polariton modes to be weU separated, apply 
another rotating wave approximation where we neglect the term Hcc and convert the 
nonlinearity term HnUn into Kerr form and get a renormalization of the polariton 
frequency and a density-density coupling between the polariton modes. This requires 
the difference in frequencies for the unperturbed modes uj^_^—uj^_ = Acj^ + involved 
to exceed the magnitude of the coupling between the modes and the nonlinearity, 

X/ACJ2 + > Jo 

Plugging in realistic values for the parameters, extracted for example from [30], {Ec = 
0.4Ghz, ijjr = 43, 6Ghz Ej/Ec = 0..150) we realize that the second inequality is indeed 
fulfilled. Engineering the capacitance Cj such that Jq is of the order of Ec/12^ we 
can ensure that the first equality is fulfilled as well. The rotating wave approximation 
eliminates the intermode exchange coupling and we obtain a Bose-Hubbard Hamiltonian 
for both modes, c+ and c_, with a density-density coupling between them, 

H^''> = H,^+H,_+Had, (10) 



where 



Hc+ = {^+4,+Ci,+ - J+ (c{+Ci+i,+ + h.c. 



^ t t 

9 



He. 



with 



_^ t t 
Hdd= -Y 2f^+-4+^i-+4-^i.- ' 

i 

UJ+ =uj'^ + Ec {cos\e) + sin2(^) cos^{e)) 
CO- =uj'_ + Ec {sm\e) + sin^{e) cos^{9)) 
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Figure 2. Plot of the effective Bose-Hubbard parameters for the c+ polariton mode, 
/7+ and J+, and the c_ polariton mode, U- and J_. As the character of the polariton 
modes change as a function of Ej / Ec from resonator field mode to qubit excitation for 
the c+ polariton modes and vice versa for the c_ polariton modes the Bose-Hubbard 
parameters change accordingly. 



J_ = Josin2(^) 
[/_ =Eccos^{e) 

We thus arrived at a two component Bose-Hubbard model for the modes c+ and c_ 
with attractive interactions and a density-density couphng between both species. The 
two species are a mixture of stripline resonator field mode and qubit excitations (9a-b) 
with different weights of the photonic or qubit contribution depending on the value of 
Ej/Ec. 

For small values of Ej/Ec the c+ polaritons become increasingly photonic. 
Consequently, their tunneling rate J+ approaches the tunneling rate of bare photons, 
Jo, and their on-site interaction vanishes. For large Ej/Ec^ on the other hand, J+ 
vanishes and the nonlinearity [/+ approaches the nonlinearity of the qubit s, Ec- For the 
c_ polaritons, the roles of both limits are interchanged. 

For each value of Ej/Ec^ the separation between the resonance frequencies of c+ 
and c_ polaritons, |cl;+ — cl;_|, is sufficiently large such that, in a scenario where we drive 
the first cavity by a microwave source, we can always adjust the frequency of the drive 
to only selectively excite one of the modes. For reasons that will become clear later we 
therefore choose the c+-polaritons to be our quantum simulator for a driven dissipative 
Bose-Hubbard model. 
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2. 5. Validity of the approximations 

To further illustrate the validity of our approximations we compare the eigenenergies 
of the full Hamiltonian i/, c.f. ([T]), approximated under assumption ([5]), with the 



eigenenergies of the Bose-Hubbard Hamiltonian H^^\ c.f. (10). The single-site 
Hamiltonians summed up in the full Hamiltonian describe the interaction between 
transmon qubit and stripline resonator in a rotating wave approximation. This 
Hamiltonian has already been used to describe an experiment revealing the nonlinear 
response of a resonator and transmon qubit system with excellent agreement between 
theory and experimental data [33] . Therefore comparison of the eigenvalues of our Bose- 
Hubbard Hamiltonian and the eigenvalues of the full Hamiltonian provides a good means 
to estimate the effects of the approximations we made. For simplicity we restricted our 
model to two sites. 

Both Hamiltonians conserve the total number of excitations and we can diagonalize 
them in each subspace with a fixed number of excitations independently. Eigenvalues 
of the full Hamiltonian in the one excitation subspace are plotted in solid lines in 
figure [3^. Without qubits, the Hamiltonian of the two resonators has eigenmodes 
a± = (ai ± a^j\f2. In figure [3] a) we also plotted the energies of these eigenmodes 
of the two coupled empty stripline resonators marked by two horizontal dash-dotted 
gray lines and the eigenenergy of the transmon qubit marked by a dash-dotted gray 
line. 

The two eigenenergies approximated by the c_-polariton mode (blue lines in figure [3] 
a) evolve from the transmon qubit energy for small values of EjjEc to the energies of 
the two stripline resonator states for large values of Ej/Ec thereby confirming our 
earlier comment that the c_-polaritons evolve from pure qubit excitations to photonic 
excitations. The eigenenergies of the c_ polaritons are degenerate for small values of 
Ej j Ec because the transmon qubits decouple from the stripline resonators and thereby 
also from each other. In figure [3] b) we plot diflFerences between the eigenenergies of the 
full Hamiltonian and the respective eigenenergies of the Bose-Hubbard Hamiltonian. For 
small values of Ej/Ec we find aberrations due to the error we make in approximating 
the Hamiltonian of the transmon qubit Hi^site.i ^i-sife,r There are also aberrations 
in the anticrossing area which are due to the neglected interactions between the c+- and 
c_-polaritons. 

For the two eigenenergies of the c+-polaritons (red lines in figure [3] a) there is a 
rather similar scenario. They approximate the two eigenenergies of the full Hamiltonian 
that are purely photonic for small Ej/Ec and evolve into qubit excitations as the 
ratio of Ej/Ec increases. There are aberrations in the anticrossing area between the 
eigenenergies of the full Hamiltonian and the Bose-Hubbard Hamiltonian for the c+- 
polaritons because of the neglected interactions between the c_- and c+-polaritons, 
plotted in figure [3] b) but there is no aberration for small values of Ej/Ec caused 
by errors made in the transmon qubit Hamiltonian because the c+-polaritons are pure 
photonic excitations for small values of Ej/ Ec- 
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Additionally to the differences between the eigenenergies of the full Hamiltonian and 
the Bose-Hubbard Hamiltonian in the one excitation subspace we plotted the differences 
in the two excitation subspace in figure |3]c). These eigenenergies can be grouped for 
the Bose-Hubbard Hamiltonian according to the distribution of excitations among the 
two polariton species. Differences of eigenenergies for states with two c_ polaritons 
are plotted in blue, for two c+ polaritons in red and for one c_ polariton and one c+ 
polariton in green. In the two excitation subspace we have similar findings as in the 
single excitation subspace. There are aberrations for the anti-crossing area because of 
the neglected intermode polariton exchange interaction. In addition, states containing 
c_ polaritons have aberrations for small values of Ej / Ec due to the approximations of 
the transmon Hamiltonian whereas c+ polaritons do not. 

Therefore the Bose-Hubbard Hamiltonian for the c+ polaritons mimics the 
behaviour of the full Hamiltonian for the full range of Ej/Ec^ provided the intersite 
coupling Jo is at most of the order of the on-site nonlinearity Ec and the polariton 
densities are not to high. To conclude: In a driven dissipative setup where we 
selectively excite the c+-Polaritons we do have a quantum simulator for a Bose-Hubbard 
Hamiltonian. 

3. Polariton statistics in the driven dissipative regime 

In this section we make use of the above explained mapping of the full Hamiltonian H 
to a two component Bose-Hubbard Hamiltonian H^^^ and consider a chain of coupled 
resonators, where we coherently drive the first resonator and adjust the microwave 
drive frequency to selectively excite the c+-polaritons. In the driven dissipative regime 
we expect to explore new physics that go beyond the equilibrium features that are 
commonly examined in many body physics. We thus calculate the polariton density 
and the density-density correlations g^'^^ in a master equation approach and analyse the 
dependencies on the system parameters J+,t/+ and the Rabi frequency of the microwave 
drive Vt. 

First experimental realisations of coupled stripline resonators are expected to 
consist of only a few resonators. To closely approximate the expected experiments and to 
speed up numerical calculations, we thus focus on a minimal chain of only two resonators. 
More specifically, we consider two stripline resonators coupled to transmon qubits, where 
the first stripline resonator is driven by a microwave source and the output signal of the 
second cavity is monitored as a function of the microwave drive frequency and the ratio of 
Ej / Ec which can be controlled by applying an external magnetic flux to the transmon 
qubits c.f. figure [4j This setup and very similar setups are currently investigated in 
experiments for example [16], and the spectroscopic measurement technique proposed 
here has already been demonstrated in single site experiments for example in [H]. 

The ouput fields are linear functions of the intra-cavity field in the second resonator 
and thus show the same particle statistics. We therefore calculated the polariton density 
and the g'^^^-function for the second cavity. To do this, we use a master equation 
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Figure 3. The accuracy of our approximations. Plotted are the first four excited 
eigenvalues of the full Hamiltonian in figure a). Eigenvalues approximated by the 
c_ polaritons are plotted in blue and eigenvalues approximated by the c+ polaritons 
in red. Figures b) and c) show the differences between the eigenenergies of the full 
Hamiltonian H and the Bose-Hubbard Hamiltonian H^^\ AE^ in the one excitation 
and two excitation subspace respectively. Differences involving eigenstates containing 
c_ polaritons are plotted in blue and differences involving eigenstates containing c+ 
polaritons are plotted in red. Differences of eigenvalues of the full Hamiltonian and 
eigenvalues of the Bose-Hubbard Hamiltonian with mixed c+ and c_ parts are plotted 
in green 



Jo 



^7 



7 



Figure 4. Sketch of experimental setup consisting of two stripline resonators coupled 
to transmon qubits with spectroscopic measurement technique: the first stripline 
resonator is excited by coherent microwave source and the output field at the second 
cavity is monitored. Stripline resonators decay at rate n mainly because of the 
finite refiectivity of the coupling capacitances and the transmon qubits decay Purcell 
enhanced at rate 7 into modes not confined by the cavity. 
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approach in which each element, the striphne resonators and the transmon qubits, couple 
to separate environments with decay rates denoted k for the stripline resonators and 7 
for the transmon qubits. Absolute values can for example be extracted from [i33j where 
7 = 3.7MHz . Decay of the stripline resonator is due to the finite transparency of the 
coupling capacitors at both ends of the resonators and decay rates for example in [30] are 
= 5.7MHz. Both environments, for the transmon qubit and the stripline resonator, are 
assumed to be in a vaccum state which is a valid assumption at typical temperatures 
for circuit QED experiments of T = 15m K. Therefore in a master equation for a 
Hamiltonian expressed in the operators for the resonator field mode a and the transmon 
qubit b the dissipators read. 



'^V[a,a]p+^V[b,b]p, 



with 



V[A, B]p = 2ApB^ - {A^Bp + pA^B) . 
These can be cast into dissipators expressed in the polariton modes c+ and c_ 



Tc+V[c+, c+]p + Tc_V[c-,C-]p + 
+ A(P[c+,c_]+P[c_,c+])p, 



with 



Kcos^{0) +7sin^(6') 
2 

Ksin^{6) +7008^(6') 



A = sin(^)cos(^) 



K — 7 



In the driven dissipative case where we selectively excite the polariton c+ mode we can 
neglect the dissipators of the c_ polaritons V[c-,C-] and the mixed dissipators T>[cj^, c_] 
and T>[c_, c+]. With these assumptions the master equation for a two site chain of the 
polariton c+-mode reads, 



dp 
di 



H. 



c+ 



f]cos(^) 



+ (^'[c+,1, C+,i]p + V[c+^2, C+,2]p) 



fill 



-,1 



+ + H, 



To solve this master equation numerically we use it to derive the coupled equations 
of motion for the expectation values of normally ordered moments of the creation and 
annihilation operators ^+,1? <^+,2 
d 



dt 



,l^+,2 



tr 



fn k J 



l^+,2 



We truncate this set of coupled equations by omitting couplings to mean values with 
n + m + k + I bigger than some rimax and solve the reduced set of equations of motion. 
To confirm the accuracy of our approach, we test its convergence with increasing nmax- 



That is, we repeat the procedure for 



+ 1, compare the results and increase 
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Ej/Ec 



Figure 5. Logarithmic density plot of the polariton density \C2C2I in the last cavity 
plotted against Ej/Ec and the frequency of the microwave drive in units of the 
frequency of the stripline resonator uo^w/ujr- Resonances in the density of polaritons 
arise where the microwave frequency matches one of the transition frequencies of the 
non-driven conservative system Hamiltonian Hc^ 



the value for n^ax in case both results differ by more than some required threshold 
value. The advantage with respect to a method that truncates the Hilbert space at 
some maximal number of excitations, is that our method becomes exact in the limit 
where the Hamiltonian becomes harmonic which is the case for small values of EjjEc- 
Moreover we experience a substantial decrease in cpu-time for this method. 

3.1. Polariton density 

We are interested in the field particle statistics in the driven dissipative regime and its 
dependencies of the on-site nonlinearity the intersite coupling J+ and the strength 
of the microwave drive Vt. We therefore first consider the density of c_ polaritons in the 
last resonator. Figure shows the density of c+ polaritons, (c\c2j^ in the second cavity 
as a function of the ratio Ej/Ec and the microwave drive frequency cj^^. The density 
of polaritons in the last cavity exhibits resonances when the microwave drive frequency 
matches one of the transition energies of the undriven conservative system Hamiltonian 
i7c+ and decreases rapidly because of the small decay rate F. One can clearly see the 
resonances due to transitions driven between the groundstate and eigenenergies in the 
one excitation subspace plotted in figure |3] a) . Transitions from the groundstate into a 
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Ej/Ec 

Figure 6. Density plot of the ^^^^ -function of the polariton mode in the second 
cavity plotted against Ej/Eq and the frequency of the microwave drive in units of 
the frequency of the stripline resonator uJ^w/^r- For increasing nonlinearity the g^'^^- 
function shows bunching regions for red-detuned microwave drive with respect to the 
energy of 1-excitation states and anti-bunching regions for blue-detuned microwave 
drive. 



two excitation state are much weaker owing to the finite Rabi frequency of the microwave 
drive ^. 

3.2. Density- density correlations 

We now consider the density-density correlations g^^^ in the last resonator. The g^^^- 
function is a quantity that describes the likelihood to measure two photons at the same 
place. The g^^^ function of the last resonator is the normalized meanvalue of the second 
order moment of the field operators in the last resonator, 



t 



Classical, thermal fields have ^'^^^-values larger or equal to unity with the coherent field 
exhibiting a ^'^^^-value of 1. A ^'^^^-value below 1, meaning that the photons are anti- 
bunched, is a sufficient condition to call the field quantum mechanical in the sense that 
there is no classical field showing the same results in measurements of the g'^^^-function. 
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With recently developed refinements of microwave measurement techniques [311 [32] . 
measurements of ^'^^^-functions in circuit QED are now becoming feasible. 

In figure [g] we plotted the ^'^^^-function of the field in the last stripline resonator. To 
get a more detailed insight of the processes leading to a ^''^^^-value for specific parameters 
we plotted the ^'^^^-function along special values of the microwave drive frequency and 
the ratio Ej / Ec marked by white lines in figure |6| Figures [s] and [9] show the results 
for the different paths, denoted by a), b), and c) in the density plot of the ^'^^^-function 
in figure ^ for the ^'^^^-function as well as the corresponding values for the density of 
polaritons in the last cavity, (c^ 2^+,2), and the second order moment, (c^ 2^+,2<^+,2C+,2)- 

For small values of Ej/Ec our system is basically linear because the nonlinearity 
= {Ec/2) sin^(0) is negligible. A harmonic field mode driven by a coherent source 
is in a coherent state. Therefore the ^'^^^-function is equal to one for small values of 
Ej/Ec' As the nonlinearity grows for increasing values of Ej/Eq the ^'^^^-function 
plotted against the ratio of Ej/Ec and the frequency of the microwave drive cj^^ gets 
more structured. In the density plot of the ^'^^^ -function in figure [g] we can identify 
resonances where the frequency of the microwave drive matches the eigenenergies of 
the unperturbed system without microwave drive and dissipation. These resonances 
manifest themselves as separating lines between bunching regions (values of g^^^ > 1) 
and anti-bunching regions (values of ^'^^^ < 1). 

To understand the origin of these separating lines, it is illustrating to analyze our 
system in terms of a symmetric mode, and an antisymmetric mode, where 

rather than the two localized modes c+^i and In terms of dj^ and c/_ the Hamiltonian 
reads, 

Hc^ = (cj+ — J+) d^_^d^ + {(jj^ + J+) d^_d- — (^d^_^d^_^d^d^ + 

- U+dld+dld_ - ^ (^dldld_d_ + /i.e.) . (12) 

The Hilbert space of the Hamiltonian i7c+ can be described by two different bases, states 
that are labeled by the number of excitations in the collective modes, 

or states that are labeled by the number of excitations in the localized modes, 

{c\ ^Y (c^ 2)"^ I \ 

|00) \nm)^ . 



The lines separating bunching and anti-bunching regions in figure [6] can now be identified 
with the energies of the 1 excitation states, 

d+|00) = |10)^^ = -L(|10)^ + |10)J (13a) 
d_|00) = |10)^^ = i=(|10)^-|10)J (136) 
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Figure 7. Sketch of the energy spectrum of the Bose-Hubbard Hamiltonian Hc_^ for 
a two site model for vanishing nonhnearity t/+ compare a) and vanishing intersite 
coupling J+ compare b). For vanishing nonhnearity a microwave drive can drive 
multiple transitions leading to a coherent state. Contrary to the linear case for strong 
nonhnearity one can only drive a transition between two distinct states as the energy 
differences between the eigenenergies aren't degenerate any more 



and the energy of a 2-excitation state, 

dldl |00) = |11)_ = (|20), - |02) J . (13c) 

To understand the origin of the anti-bunching regions for a microwave drive that is blue 
detuned with respect to the energies of the states ( 13a[ -c) and the bunching regions for a 
red detuned microwave drive one has to consider the spectrum of the Hamiltonian . 

For small nonhnearity, that is for values of Ej/Ec < 50, the Hamiltonian (12) 
reduces to a Hamiltonian for two uncoupled harmonic oscillators described by the modes 
dj^ and with energies ujj^ — J+ and ujj^ + J+ respectively. The eigenenergies in this 
situation are shown in figure [7] a). A microwave drive with frequency ujj^ — J+ as depicted 
in figure [7| not only drives the transition from the groundstate to the first excited state of 
the symmetric collective mode |0 0)^^ |1 0)cm ^^^^ other transitions to higher 
excited states \nQ)^^ \n-\- 10)^^. As a result the steady state in this situation is 
always the coherent state exhibiting a ^'^^^-value of 1. For slightly increased values of 
the nonhnearity that remain in the range [/+ < Fc^, the system can still be described 
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in terms of two weakly interacting collective modes. But the symmetric as well as the 
antisymmetric mode are subject to the nonlinearity and an intermode interaction, c.f. 



Hamiltonian (12). This can be seen in figure [s] a) where we plotted the ^'^^^-values 
that deviate from the value of a coherent field. The ^'^^^ -function shows anti-bunching 
regions for blue detuned microwave drive with respect to the energies of the states 
(13a-c) and bunching regions for red detuned microwave drive. To gain insight into 
the underlying physical principles in this situation we calculated the density (c^2^+,2) 
and the second symmetric moment (c^ 2^+ 2^+,2C+,2) by an iterative meanfield approach 
to solve the master equation (11) with the Hamiltonian written as in (12). Operator 
mean values of a single driven dissipative mode with Kerr nonlinearity can be computed 
exactly |[34j and we expand this model in a meanfield way to incorporate the denisty- 
density coupling. With this method we get good agreement with the numerical exact 
values for the density in the last cavity and are able to compute values for the polariton 
density close to the systems eigenenergies ( 13a[ -c) where our numerical approach fails to 
converge. For details about the method please see [Appendix A[ These results support 
our assertion that the system can be described by weakly interacting collective modes in 
the limit of small nonlinearities In figure [s] a) numerically exact values are plotted 
in solid lines and values obtained by the above mentioned mean field method are plotted 
in dashed lines. 

For strong nonlinearity and small intersite coupling that is for values of 
Ej/Ec > 50, the c+-polaritons become transmon excitations and the Hamiltonian 
splits into two parts describing the first and the second transmon qubit respectively. 
Here the collective modes and (i_ no longer decouple and the localized modes Ci and 
C2 become a more appropriate description of the system. The eigenenergy spectrum 
in this situation is shown in figure [T] b) . The main diflFerence to the spectrum without 
nonlinearity is that the microwave drive can't be adjusted to drive multiple transitions. 
In order to drive the transition to the state |02)^ for example one has to adjust the 
microwave frequency to match half of the energy difference between the groundstate 
and the 2-excitation state |02)^ because it is a two photon transition. Due to the 
anharmonicity of the eigenenergy spectrum no other transition can be driven. The 
difference of microwave frequencies needed to drive the transition from groundstate to 
|01)^ respectively |02)^ amounts to which is bigger than the linewidth Fc^. To 

get an estimate for the value of g^'^'^ we simplify our model assuming that the frequency 
of the microwave drive is adjusted such that it resonantly drives a transition between 
the groundstate of our model \00)^ and some excited state |0n)^. Provided the Rabi 
frequency Q and the loss rates and 7 are all small compared to the frequency separation 
between different resonance lines, the system can then be modeled by a two level system 
consisting of the groundstate of our model |00)^ and the excited state |0n)^. In this 
situation the maximal occupation inversion one could get in the steady state is, 

Pmax = 1(100) (00| + |0n) (On|) , 
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and the ^r'-^-'-value for this density matrix would be 

2n(n - 1) 



tr 



tr 



which is below one for a 1-excitation state and above one for every state containing 
more than 2 excitations. Therefore bunching areas arise if states containing more than 
2-exciations are excited and anti-bunching areas arise if only 1-exitation states can be 
excited and the photons pass the setup "one by one" . In our Bose-Hubbard model the 
on-site nonlinearity is negative and hence all transition frequencies to states containing 
more than two excitations are red detuned with respect to transition frequencies to states 
containing only one excitation ( 13a| -c). This is why bunching areas arise for red detuned 
microwave drive and anti bunching areas arise for blue detuned microwave drive. 

If we adjust the microwave drive frequency for every value of Ej/Ec to match 
the eigenfrequency of the antisymmetric 1-excitation state we get the transition from 
a perfectly uncorrelated field with g^'^^ = 1 to strongly correlated, anti-bunched field 
statistics with g^'^^ < 1 see figure [o] For a quantum phase transition of the ground state 
of the Bose-Hubbard Hamiltonian one would expect this transition as a consequence 
of the interplay of the intersite hopping J+ and the on-site nonlinearity For the 
driven dissipative system we observe that the interplay between the Rabi frequency of 
the microwave drive Qco8{9) and the on-site nonlinearity determines the particle 
statistics. This can be seen in figure [o] where we plotted the ^'^^^-function and the 
intersite coupling, on-site nonlinearity and the Rabi frequency of the microwave drive. 



4. Summary 

We have shown that a chain of capacitively coupled stripline resonators each coupled to 
a transmon qubit can be described by a Bose-Hubbard Hamiltonian for two species of 
polaritons. The validity of our approach has been checked for realistic parameters of the 
transmon qubits and stripline resonators and for low densities of polaritons. In a driven 
dissipative regime where a microwave source coherently drives the first cavity one can 
selectively excite only one species of polaritons and investigate the properties of a driven 
dissipative Bose-Hubbard model. We calculated the density and the ^'^^^-function of the 
polaritons in the last resonator of a two site setup and investigated their dependencies on 
the microwave drive, the intersite coupling and the on-site nonlinearity. For vanishing 
nonlinearity the ^'^^^-function is approximately equal to unity indicating a coherent 
field. With increasing nonlinearity bunching and anti-bunching areas arise depending 
on the frequency of the microwave drive. For a microwave drive that is in resonance 
with a transition to a state with one excitation, the polaritons are anti-bunched. If, 
on the other hand, the microwave drive can resonantly excite states containing more 
than 2 excitations in a multi-photon transition, the polaritons become bunched. If we 
adjust the microwave drive frequency to match one of the system's single excitation 
eigenenergies and compute the ^'^^^ -function for different values of on-site nonlinearity. 



CONTENTS 



20 



1.1 

1.05 
1 

0.95 



0.9 
0.95 



a) Ej/Ec = 25 




J 










^ " 



1.05 




0.95 



1 1.05 



1.1 



1.27 



b) Ej/Ec = 125 



1.28 



1.29 



1.3 



10^ 








10-^ 






(4^2) 


10-^ 

1 1 
10° 








27 


1.28 1.29 


1 


3 










10-^ 




1 ' A' 






- - icW^c,)^^ 




1 In 1 / II 






10-" 








\ 1-ex. transition 


10-^ 


^111 X/l 1 1 1 




j 2-ex. transition 


10-^ 


III II II 




j 3-ex. transition 



Figure 8. Plots of the ^^^^-function, the density of polaritons \c\c2i and the second 

order moment (^c\c2C2C2^ in the second cavity for special values of Ej /Ec are shown. 
For all plots we have chosen the intersite coupling constant and the on-site nonlinearity 
to be J^/uOr = 0.04 and Ec/oOr = 0.04 and the decay rates of transmon qubit and 
stripline resonator to be 7/6^^ = 0.00008 and n/uor = 0.00004 respectively but we 
applied different Rabi frequencies of the microwave drive: for a) Vt/uOr = 0.004, and 
for b) Vt/uJr = 0.001. Plotted are results obtained by numerical calculation of the 
masterequation in solid lines and results obtained by a mean field approach with 
an exact single-site solution in dashed lines. Eigenenergies of the system without 
dissipation and driving are signalized by vertical dash-dotted lines. For Ej/Ec = 25 
one can see clearly separated resonances for the symmetric and antisymmetric states 
d}j^ 1 00) and a two photon resonance for the state d}j^d^_ |00). The shape of the 
resonances at d^j^ |00) and d^_ |00) are reproduced by the meanfield approximation and 
therefor a single mode feature. The two photon resonance for the sate d^'^dL |00) 
is not reproduced by the meanfield approach since it does not correctly incorporate 
the interactions between (i^ and (iL modes. For Ej/Ec = 125 multiple resonances 
determined by the eigenenergies of the system without dissipation and driving arise . 



intersite coupling and Rabi frequency of the microwave drive we see a transition from 
coherent to anti-bunched field statistics. That is, the polaritons are uncorrelated for 
small nonlinearity and exhibit a transition to anti-bunched behaviour as the on-site 
nonlinearity becomes larger than the Rabi frequency of the microwave drive. All our 
findings could be explored in experiments based on readily available technology. 
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Figure 9. A plot of the ^^^^-function, the density of polaritons \ c\c2i and the second 

order moment ^C2C2C2C2^ in the second cavity is shown. The microwave drive frequency 
is chosen to match the transition from the ground state to the symmetric 1-excitation 
eigenstate. For ah plots we have chosen the intersite coupling constant and the on-site 
nonlinearity to be J^/uOr = 0.04 and Ecl^r = 0.04, the decay rates for transmon qubit 
and stripline resonator to be 7/co'r = 0.00004 and njujr = 0.00008 respectively and 
the Rabi frequency of the microwave drive to be VtjuJr = T/ujr = 0.00004. Plotted 
are results obtained by numerical calculation of the masterequation in solid lines and 
results obtained by a mean field approach with an exact single-site solution in dashed 
lines. In resonance to the symmetric state ^^^^ shows a transition from uncorrelated 
coherent field particle statistics to anti-bunched correlated field particle statistics. The 
transition from coherent to anti-bunched is determined by the interplay of the Rabi 
frequency of the microwave drive and the on-site nonlinearity. 
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Appendix A. Meanfield approximation for two coupled modes with small 
nonlinearity 



We want to solve the master equation (11) for smaU values of Ej/Ec- For Ej/Eq = 25, 
we find J+ ^ and express Hamiltonian i7c+ iii terms of the collective modes and 



ex 
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where, 

Ha^ = ^1^^ (^dl + + (Ac^^ ± j^) ^^^^ _ ^dUid±d^ 

Hdd = — tZ+d^rf+rf^d- 

Hex = - ^ {d\d\d_d_ + /i.c.^ . 

Here Acl;+ = cl;+ — cj/ is the difference between the frequency ujj^ and the microwave 
drive frequency uji. Since J+ the two coUective modes are energeticaUy far 

separated and we neglect H^x in a rotating wave approximation. In zeroth order in our 
approximation we neglect Hdd as well and decouple the collective modes completely. The 
driven dissipative master equation for the collective modes rf+ and rf_ with Hamiltonians 
Hd^ and Hd_ can be solved exactly for steady state values, c.f. [34j. We compute the 
densities of the two collective modes = in zeroth order approximation and 

use them to approximate the intermode density-density coupling, 

2 

This way the modes are still decoupled but the density of the rf+-mode induces a 
frequency shift of the rf_-mode cl;+ — J+ ^ cj+ — J+ — and vice versa. We proceed 

to calculate the densities of the two modes with the shifted frequencies to obtain the 
densities of the two modes in first order approximation We iterate this method 
until the diflFerence between densities of consecutive order in approximation is smaller 
than some threshold value. After this procedure we approximate the density in the last 
stripline resonator 

4,2C+,2) ^ \ {{dU+) + (did.)) + Re [(4) (d.) 

where we calculated the values for and with the renormalized frequencies 

cl;+ + J+ — t/^^f^^^^) and uj^ — — [7+5^'^''"'^ With this procedure we can reproduce 
the numerically exact values for the density in the last stripline resonator and are able 
to compute values for the polariton density close to the systems eigenenergies 
where the numerically exact method converges very slowly and becomes numerically 
very demanding. 
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